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Abstract. 

A possible candidate for the present accelerated expansion of the Universe is "phantom energy", which possesses an 
equation of state of the form £0 = p/p<— 1, consequently violating the null energy condition. As this is the fundamental 
ingredient to sustain traversable wormholes, this cosmic fluid presents us with a natural scenario for the existence of these 
exotic geometries. In this context, we shall construct phantom wormhole geometries by matching an interior wormhole 
solution, governed by the phantom energy equation of state, to an exterior vacuum at a junction interface. Several physical 
04 ' properties and characteristics of these solutions are further investigated. The dynamical stability of the transition layer of 

these phantom wormholes to linearized spherically symmetric radial perturbations about static equilibrium solutions is also 
explored. It is found that the respective stable equilibrium configurations may be increased by strategically varying the 
wormhole throat radius. 
> ■ 

O 1. INTRODUCTION 

cn ^ 

, One of the most challenging current problems in cosmology is explaining the accelerated phase of expansion of the 

' Universe [1]. Several candidates, responsible for this cosmic expansion, have been proposed in the literature, namely, 

\ dark energy models, generalizations of the Chaplygin gas, modified gravity and scalar-tensor theories, tachyon scalar 

O ■ fields and braneworld models, such as the Dvali-Gabadadze-PoiTati (DGP) model, amongst others. The dark energy 

models are parametrized by an equation of state given by © = p/p, where p is the spatially homogeneous pressure 
5-H . and p is the dark energy density. For the cosmic expansion, a value of O < — 1 /3 is required, as dictated by the 

Friedman equation a/ a = —4n{p + p/3), so that a > 0. A specific exotic form of dark energy, denoted phantom 
energy, has also been proposed, possessing the peculiar property of O) < — 1 [2] . This parameter range is not excluded 

■ by observation, and possesses peculiar properties, such as the violation of the null energy condition and an infinitely 
, increasing energy density, resulting in a Big Rip, at which point the Universe blows up in a finite time [2]. However, 

■ recent fits to supernovae, CMB and weak gravitational lensing data indicate that an evolving equation of state co 
crossing the phantom divide —1, is mildly favored [3]. 

As the phantom energy equation of state violates the null energy condition, p + p < 0, which is the fundamental 
ingredient to sustain traversable wormhole [4, 5], one now has at hand a possible source for these exotic spacetimes. 
However, a subtlety needs to be pointed out, as emphasized in Refs. [6, 7]. The notion of phantom energy is that of 
a homogeneously distributed fluid. When extended to inhomogeneous spherically symmetric spacetimes, the pressure 
appearing in the equation of state is now a radial pressure, and the transverse pressure is then determined via the field 
equations. Wormhole geometries were found in Ref. [6], by considering specific choices for the distribution of the 
energy density, and in Ref. [7], a complementary approach was traced out, by imposing appropriate choices for the 
form function and/or the redshift function, and the stress-energy tensor components were consequently determined. 
This latter approach shall be traced out in this paper. 

It is also of a fundamental importance to investigate the stability of these phantom wormhole geometries. We 
shall model these spacetimes by matching an interior traversable wormhole geometry with an exterior Schwarzschild 
vacuum solution at a junction interface [8]. In this work, we analyze the stability of these phantom wormholes to 
linearized perturbations around static solutions. Work along these lines was done by considering thin-shell wormholes, 
using the cut-and-paste technique [9]. The advantage of this analysis resides in using a parametrization of the stability 
of equilibrium, so that there is no need to specify a surface equation of state. The stability analysis of these thin-shell 
wormholes to linearized spherically symmetric perturbations about static equilibrium solutions was carried out by 
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assuming that the shells remain transparent under perturbation [lOJ. This amounts to considering specific spacetimes 
that do not contribute with the momentum flux term in the conservation identity, which provides the conservation law 
for the surface stress-energy tensor. The inclusion of this term, corresponding to the discontinuity of the momentum 
impinging on the shell, severely complicates the analysis. However, we shall follow the approach of Ishak and Lake 
[10], with the respective inclusion of the momentum flux term and find stability equilibrium configurations of specific 
phantom wormhole geometries. 

This paper is organized in the following manner. In Section 2, we outline the general properties of phantom energy 
traversable wormhole, and in Section 3, we present specific solutions of these phantom wormholes. We also show 
that using these specific constructions, and taking into account the "volume integral quantifier", one may theoretically 
construct these spacetimes with infinitesimal amounts of ANEC violating phantom energy. In Section 4, we outline 
a general linearized stability analysis procedure, and then apply this analysis to phantom wormhole geometries and 
determine their respective stability regions. Finally in Section 5, we conclude. 



2. PHANTOM ENERGY TRAVERSABLE WORMHOLES 

The spacetime metric representing a spherically symmetric and static wormhole is given by 

ds" = dt^ + —^frrr + {dd^ + sin' Od(^^), (1) 

1 —b{r)/r 

where 4>(r) and b{r) are arbitrary functions of the radial coordinate, r. 4>(r) is denoted as the redshift function, for 
it is related to the gravitational redshift; b{r) is called the form function [4]. The radial coordinate has a range that 
increases from a minimum value at tq, corresponding to the wormhole throat, to a, where the interior spacetime will 
be joined to an exterior vacuum solution. 

Using the Einstein field equation, G^y = Stt T^y, in an orthonormal reference frame, (with c = G = 1) we obtain the 
following stress-energy scenario 

b' = Snr^p, (2) 
b + Snr^Pr 

2r^{l-b/r)' ^' 

p'r = ^{Pt-Pr)-{P + Pr)^', (4) 

where the prime denotes a derivative with respect to the radial coordinate, r. p(r) is the energy density, pr{r) is the 
radial pressure, and pt{r) is the lateral pressure measured in the orthogonal direction to the radial direction. Equation 
(4) may be obtained using the conservation of the stress-energy tensor, J^.y = 0. At the throat we have the flaring out 
condition given by {b — b'r)/b^ > [4, 5], which may be deduced from the mathematics of embedding. From this 
we verify that at the throat b{ro) = r = ro, the condition b'iro) < 1 is imposed to have wormhole solutions. For the 
wormhole to be traversable, one must demand that there are no horizons present, which are identified as the surfaces 
with e^* —^ 0, so that <J>(r) must be finite everywhere. Note that the condition 1 — fo/r>Ois also imposed. 

A fundamental ingredient of traversable wormholes and phantom energy is the violation of the null energy condition 
(NEC). Matter that violates the NEC is denoted as exotic matter. Note that the notion of phantom energy is that of a 
homogeneously distributed fluid; however, when extended to inhomogeneous spherically symmetric spacetimes, the 
pressure appearing in the equation of state is now a radial pressure, and the transverse pressure is then determined 
from Eq. (4). Using the equation of state representing phantom energy, pr — (op with co<—l, and taking into account 
Eqs. (2)-(3), we have the following condition 

, b + corb' 
^^'^ = 2rHl-b/r) - 

We now have four equations, namely the field equations, i.e., Eqs. (2)-(4), and Eq. (5), with five unknown functions of 
r, i.e., p{r), Pr{r), Pt{r), b{r) and ^(r). To construct solutions with the properties and characteristics of wormholes, 
we consider restricted choices for b{r) and/or 4>(r). In cosmology the energy density related to the phantom energy 
is considered positive, p > 0, so we shall maintain this condition. This impUes that only form functions of the type 
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b'{r) > are allowed. We can construct asymptotically flat spacetimes, in which b{r)/r ^ and <I> — > as r — > oo. 
However, one may also construct solutions with a cut-off of the stress-energy, by matching the interior solution to an 
exterior vacuum spacetime, at a junction interface. 



3. SPECIFIC PHANTOM WORMHOLE MODELS 
3.1. Asymptotically flat spacetimes 

To construct an asymptotically flat wormhole solution |7|, consider = <I>o = const. Thus, from Eq. (5) one 
obtains 

b{r)=n,{r/r^)-"'', (6) 

so that fe(r)/r = (ro/r)(i+®)/® ^ for r ^ oo. We also verify that b'{r) = -(1 /aj)(r/ro)"(^+'")/'*', so that at the throat 
the condition fo'(ro) = l/|ci)| < 1 is satisfied. 

The stress-energy tensor components are given by 

MO = a>pW = ^^P)"^ and MO = ^(^)p)"V (7) 
STTrQ V r / loTTrQ \ CO J \ r J 

Thus, determining the parameter oj from observational cosmology, assuming the existence of phantom energy, one 
may theoretically construct traversable phantom wormholes by considering the above-mentioned form function and a 
constant redshift function. 

Using the "volume integral quantifier", which provides information about the "total amount" of averaged null energy 
condition (ANEC) violating matter in the spacetime (see Ref . [1 1] for details), given by /y = / [p (r) -|- Pr{r)] dV, with 
a cut-off of the stress-energy at a, we have 



Iv= Hr-b) 

Jrn 



In 



.2* 



l-b/r 



dr. (8) 



Considering the specific choices for the form function and redshift function for the traversable wormhole given 
above, with o) = — 2 for simpUcity, the volume integral assumes the value ly = ro(l — ^Ja|r^. Taking the limit a ro, 
one verifies that ly = Jip + Pr) dV 0. Thus, as in the examples presented in [11], one may construct a traversable 
wormhole with arbitrarily small quantities of ANEC violating phantom energy. Although this result is not unexpected 
it is certainly a fascinating prospect that an advanced civilization may probably construct and sustain a wormhole with 
vanishing amounts of the material that comprises of approximately 70% of the constitution of the Universe. 



3.2. Isotropic pressure 



Consider now isotropic pressures, pr = pt = p- By taking into account the form function, considered above, given 
byZ7(r) = ro(r/ro)-i/'',and using Eq. (4) and Eq. (2), one finds that the redshift function is given by 



<J.(r) 



3(0+1 

T 



CO 



In 



(See Ref. [7] for details). The stress-energy tensor components are provided by 



Pir) 



(Op{r) = - 



1 



(t 



(9) 



(10) 



Note that the spacetime is not asymptotically flat. However, one may match the interior wormhole solution to an 

exterior vacuum spacetime at a finite junction surface. Using the "volume integral quantifier", Eq. (8), with a cut-off 
of the stress-energy at a, and taking into account the choices for the form function and redshift function considered 
above, with O) = — 2, the volume integral takes the following value 



Iv =a(^l- \/ro/a^ ]n 



10 



+ ^lOa-h llro-21ro-\/«/'"0 ) + « (■\/'"o/«- 1 ) In 



21/2 



(11) 
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Once again taking the limit ro, one verifies that ly — » 0, and as before one may construct a traversable wormhole 
with arbitrarily small quantities of ANEC violating phantom energy. 



4. STABILITY ANALYSIS 
4.1. Junction conditions 

We shall model specific phantom wormholes by matching an interior traversal wormhole geometry, satisfying the 
equation of state pr = (Op with (o < —1, with an exterior Schwarzschild solution at a junction interface E, situated 
outside the event horizon, a > ri, = 2M. 

Using the Darmois-Israel formalism [12], the surface stress-energy tensor, S'j, at the junction interface Z, provide 
us with the following expressions for the surface stresses 



1 

47ta 



1 \-d^ 

a 



b(a) 



1 



1 



+ (l+a<I>')(l-f+fl2 



.7.7 i^(b-b'a) 
2ia-b) 



(12) 



(13) 



where o and ^ are the surface energy density and the tangential surface pressure, respectively (see Refs. ([13, 14, 15]) 

for details). 

r .. 1 + 

where is the unit normal 4— vector to E, and 



We shall use the conservation identity, given by 5^1^. = 



e|^.j are the components of the holonomic basis vectors tangent to E. The momentum flux term in the right hand side 
corresponds to the net discontinuity in the momentum which impinges on the shell. Using = — [a + 2d{a + ^)/a], 
the conservation identity provides us with 



with 



Ana?- 



b'a-b 

a I 



2«(1-|) 



(14) 



Equation (14) evaluated at a static solution ai^, shall play a fundamental role in determining the stability regions. 

Equation (12) may be rearranged to provide the thin shell's equation of motion, i.e., c? + V (a) = 0, with the potential 
given by 

where = Anoc?- is the surface mass of the thin shell. For computational purposes and notational convenience, we 
define the following factors 



fo(a)/2+M , , , M -bid) 12 
F{a) = 1 - and G{a) = , 



(16) 



Linearizing around a stable solution situated at oq, we consider a Taylor expansion of V{a) around oq to second 
order, given by 



1. 



V{a) = y (ao) + y'(ao) {a - ao) + -V"{ao){a - aof + 0[{a- aof] . 



(17) 



Evaluated at the static solution, at a = ao, we verify that y(ao) = and y'(flo) = 0. From the condition y'(ao) = 0, 
one extracts the following useful equilibrium relationship 



2^0 



ao 



F' -2l ""^ ^ [ 



(18) 
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which will be used in determining the master equation, responsible for dictating the stable equilibrium configurations. 

The solution is stable if and only if V{a) has a local minimum at ao and V"{ao) > is verified. The latter condition 
takes the form 



~1 



aG 

Ms 



aG\ ( aO^ " 



(19) 



Using my = ^na^a, and taking into account the radial derivative of a' , Eq. (14) can be rearranged to provide the 
following relationship 

=T-4;r(7'Tj, with T = — {a + ^) + 2naZ' , (20) 

and the parameter r\ is defined as r; = S^' jo' . 

Equation (20) will play a fundamental role in determining the stability regions of the respective solutions. Note that 
the parameter ^Jy\ is normally interpreted as the speed of sound, so that one would expect that < t] < 1, based on the 
requirement that the speed of sound should not exceed the speed of light. However, in the presence of exotic matter 
this cannot naively de done so. Therefore, in this work the above range will be relaxed. We refer the reader to Ref. [9] 
for an extensive discussion on the respective physical interpretation of J] in the presence of exotic matter. 

We shall use 7] as a parametrization of the stable equiUbrium, so that there is no need to specify a surface equation 
of state. Thus, substituting Eq. (20) into Eq. (19), one deduces the master equation given by 

(7'm,Tjo>0, with 0=|i(r2-^) + i-m,T, (21) 

and Tjo = T] (ao). The master equation dictate stable equilibrium regions, given by the following inequalities 

Tjo > 0, if a'm, >0, (22) 
Tjo < 0, if (t'ot, <0, (23) 

with the definition = ©/(ff'mv). 

We shall now consider the phantom wormhole geometries found in Secion 3, and consequently determine the 
stabiUty regions dictated by the inequalities (22)-(23). In the specific cases that follow, the expUcit form of T, I*, 
and are extremely messy, so that as in [10], we find it more instructive to show the stability regions graphically. 



4.2. Stability regions 

4.2.1. Asymptotically flat spacetimes 

Consider the specific choices for the redshift and form functions, for an asymptotically flat spacetime, given by 
<I>(r) = <I>o = const and the form function given by Eq. (6), b{r) = rQ{r/ ro)^^^'^ . To determine the stability regions 
of this solution, we shall separate the cases of b{ao) < 2M and b{ao) > 2M. From Eq. (12) and the definition of 
nis — Aiia^a, this corresponds to mj > and nis < 0, respectively. Here, we shall relax the condition that the surface 
energy density be positive, as in considering traversable wormhole geometries, one is already dealing with exotic 
matter. Note that for a < 0, the weak energy condition is readily violated. 

For b{ao) < 2M, i.e., for a positive surface energy density, we need to impose the condition < 2M, so that the 
junction radius hes outside the event horizon, ao > 2M. Thus, the junction radius lies in the following range 

2M < ao < 2M(2M/ro)"('+'") . (24) 

For a fixed value of O, we verify that as ro 0, then ao °°- The range decreases, i.e., ao ~^ 2M, as rg 2M. Note 
that by fixing and decreasing O), the range of ao is also significantly increased. 

For a fixed value of the parameter, for instance (O = —2, we shall consider the following cases: ro/M = 1.0, so that 
2 < ao/M < 4; and for ro/M = 0.25, we have 2 < ao/M < 16. The respective stabihty regions are depicted in Fig. 1. 
From Eq. (14), one may show that a' < 0, and as a > 0, this implies Wja' < 0. Thus, the stability regions, dictated 
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FIGURE 1. Plots for a positive surface energy density, i.e., b(ao) < 2M. We have defined a = oq/M, and considered (O = —1 
for both cases. The first plot is given by r^/M = 1, and the second by tq/M = 0.25. The stabiUty regions are given below the solid 
curve. See the text for details. 




FIGURE 2. Plots for a negative surface energy density, considering tq/M < 2. We have defined a = ao/M, and considered 
to = — 2 for both cases. The first plot is given by ro/M = 0.5, and the second by ro/M = 1.5. The stabihty regions are given above 
the first solid curve, and below the second solid curve. See the text for details. 



by the inequality (23), lie beneath the solid lines in the plots of Fig. 1 . For decreasing values of ro/M, as oq increases, 
the values of 770 are further restricted. Using positive surface energy densities, stable geometries may be found well 
within the bound of < 770 < 1, and the stability regions increase for increasing values of ro/M. 

For b{a(j) > 2M, the surface mass of the thin shell is negative, mj (flo) < 0. We shall separate the cases of ro < 2M 

and ro > 2M. If ro < 2M, the range of the junction radius is given by gq > 2M(2M/ro) For this specific case, 

a' possesses one real positive root, R, in the above range, signalling the presence of an asymptote, a' |r = 0. We verify 
that a' < for 2M(2M/ro)"(i+®) < ao < /?, and d' > for ao > R. Thus, the stabihty regions are given by 

T/o > 0, if 2M(2M/ro)"('+'") <ao</?, (25) 
Tjo < 0, if ao>R. (26) 

Consider for ft) = -2, the particular cases of ro/M = 0.5, so that aq /M > 8, and ro/M = 1 .5, so that ao /M > 2.667. 
The asymptotes, g'\r = 0, for these cases exist at R/M ~ 13.9 and R/M o:^ 4.24, respectively. These cases are 
represented in Fig. 2. Note that for increasing values of ro /M, the range of ao decreases, and the values of 770 are 
less restricted. Thus, one may conclude that the stability regions increase, for increasing values of ro/M. 

If ro > 2M, then obviously ao > ro. We verify that a' > 0, and consequently a' < 0, so that the stabihty region is 
given by inequality (23). We verify that the values of t/o are always negative. However, by increasing ro/M, the values 
of Tjo become less restricted, and the range of ao decreases. 
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FIGURE 3. Plots for an isotropic pressure phantom wormhole. We have defined a = oq/M and considered (O = —2 for both 
cases. For b{aQ) < 2M, the condition ro < 2M is imposed. The first plot is given by ro/M = 1.0, and the second by ro/M = 1.5. 
The stabiUty regions are given below the solid curves. 



4.2.2. Isotropic pressure, Pr= Pt = P 

Consider an isotropic pressure phantom wormhole geometry, with the redshift function given by Eq. (9), •t>(r) = 
(3co + 1)/(1 + CO) ln(r/ro) and b{r) = r^ (r/ro)"'/'". To determine the stability regions of this solution, as in the 
previous case, we shall separate the cases of b{ao) < 2M and b{aQ) > 2M. 

For b{ao) < 2M, we have Ms > 0, and the condition ro < 2M is imposed. Therefore, the junction radius lies in the 
same range as the previous case, i.e., Eq. (24). We also verify that cr' < in the respective range. Thus the stability 
region is given by 

Tjo<0, if 2M<ao<2M(2M/ro)"(^+'"^ (27) 

Consider, for simplicity, CO — —2, and the cases for ro /M = 1 and ro/M = 1 .5 are analyzed in Fig. 3. The ranges are 
given by 2 < ao/M < 4 and 2 < ao/M < 2.667, respectively. Note that as ro/M decreases, the range of ao increases. 
However, the values of the parameter t]o become more restricted. Thus, one may conclude that the stability regions 
increase, as ro/M increases. 

For b{ao) > 2M, then (ao) < 0. As before, we shall separate the cases of ro < 2M and ro > 2M. For ro < 2M, the 
range of ao is given by ao > 2M(2M/ ro)"^^"'"'"^ as in the previous case of the asymptotically flat wormhole spacetime. 

For this case a' also possesses one real positive root, R, in the respective range. We have a' < for (2M/ro)^'/''' < 
a <R, and cr' > for ao > R. The stability regions are also given by the conditions (25)-(26). We have considered 
the specific cases of ro/M = 1 so that the respective range is ao/M > 4; and ro/M = 1.5, so that ao/M > 2.667. The 
asymptotes, c'Ir = 0, for these cases exist at R/M ~ 6.72 and R/M ~ 4.24, respectively. This analysis is depicted in 
the plots of Fig. 4. For decreasing values of ro/M, the stability parameter 770 becomes less restricted and the range of 
the junction radius increases. Thus, one may conclude that the stability regions increase for decreasing values of ro/M. 

5. CONCLUSION 

As the Universe is probably constituted of approximately 70% of null energy condition violating phantom energy, 
this cosmic fluid may be used as a possible source to theoretically construct traversable wormholes. In fact, it was 
found that infinitesimal amounts of phantom energy may support traversable wormholes [7]. In this paper, we have 
modelled phantom wormholes by matching an interior traversable wormhole geometry, satisfying the equation of state 
p = cop with £0 < —1, to an exterior vacuum solution at a finite junction interface. We have analyzed the stability 
of these phantom wormholes, an issue of fundamental importance, to linearized perturbations around static solutions, 
by including the momentum flux term in the conservation identity. We have considered two particularly interesting 
cases, namely, that of an asymptotically flat spacetime, and that of an isotropic pressure wormhole geometry. We have 
separated the cases of positive and negative surface energy densities and found that the stable equiUbrium regions may 
be significantly increased by strategically varying the wormhole throat. As we are considering exotic matter, we have 
relaxed the condition < tjq < 1, and found stability regions for phantom wormholes well beyond this range. 
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FIGURE 4. Plots for an isotropic pressure pliantom wormliole, for fc(ao) > 2M and rQ < 2M. We liave defined a = oq/M and 
considered a — —2 for botii cases. Tfie first plot is given by ro/M = 1.0, and the second by r^/M = 1.5. The stability regions are 
given above the first solid curve, and below the second solid curve. 



An interesting feature of the phantom regime is that due to the fact of the accelerated expansion of the Universe, 
macroscopic wormholes could naturally be grown from the submicroscopic constructions that originally pervaded 
the quantum foam [16]. The accretion of phantom energy onto Morris-Thorne wormholes was further explored in 
Refs. [17], and it was shown that this accretion gradually increases the wormhole throat which eventually overtakes 
the accelerated expansion of the universe, consequently engulfing the entire Universe, and becomes infinite at a time 
in the future before the big rip. This process was dubbed the "Big Trip" [17]. Traversable wormholes supported by 
the generalized Chaplygin gas were also found [18], and in the context of the "Big Trip", it was found that the latter 
may be avoided in a wide region of the Chaplygin parameters, by the accretion of a generalized Chaplygin gas onto 
wormholes [19]. In concluding, it is noteworthy the relative ease with which one may theoretically construct traversable 
wormholes with the exotic fluid equations of state used in cosmology to explain the present accelerated expansion of 
the Universe. 
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